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We investigate the frustrated magnetic interactions in cubic transition 
metal oxides with orbital degeneracy. The Cg orbitals order easier and 
their ordering explains the A-type antiferromagnetic phase in KCuFs and 
LaMnOa. In t2g systems the magnetic order changes at a transition from 
an orbital liquid to orbital ordered states. The fluctuations of t2g orbitals 
play a prominent role in LaVOa and YVO3, where they compete with the 
Jahn- Teller effect and trigger the C-type antiferromagnetic order. 

1. Spin-orbital physics in transition metal oxides 

Large on-site Coulomb interactions (xU m. transition metal oxides sup- 
press charge fluctuations and lead to the (partial) localization of d electrons 
which interact by the effective superexchange interactions. When such local- 
ized electrons occupy degenerate orbital states, one has to consider orbital 
degrees of freedom at equal footing with electron spins [1]. The importance 
of the orbital degrees of freedom in such systems has been emphasized long 
ago for cuprates [2] and for V2O3 [3], when it was also realized that fer- 
romagnetic (FM) superexchange could be induced by the Hund's exchange 
interaction oc Jh [4], but only recently it has been fully appreciated that 
the orbital physics leads to several novel and interesting phenomena. 

The superexchange which involves the orbital degrees of freedom is de- 
scribed by the so-called spin-orbital models [5], and is typically highly frus- 
trated even on a cubic lattice [6] . Although this frustration might even lead 
to the collapse of magnetic (or orbital) long-range order in the limit of weak 
cx Jh, in real eg systems it is largely suppressed by Jh/U ~ 0.12 [7], where 
U is the intraorbital interaction, and structural phase transitions stabilize 
a particular ordering of occupied orbitals, supporting the 74-type antiferro- 
magnetic (AF) order. Here we show that this happens even in the absence 
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of the Jahn- Teller (JT) effect in the Cg systems with degenerate orbitals 
filled either by one hole (KCuFa) [8], or by one electron (LaMnOs) [9]. 

The transition metal oxides with partly filled t2g orbitals are even more 
fascinating. The quantum phenomena are here more important and stabilize 
the coherent orbital liquid ground state in the spin 5 = 1/2 Mott-insulator 
LaTiOa [10] , which preserves the cubic symmetry and explains the observed 
isotropic G-type AF order [11]. In vanadium compounds rather involved 
spin-orbital models, which describe coexisting AF and FM interaction, were 
recently introduced for LiV02 [12] and V2O3 [13]. The superexchange is 
again frustrated in cubic systems, and C-type of AF order, observed both in 
LaVOa and T = [14] and in YVO3 77 < T < 114 K [15], can be explained 
as supported by quantum one-dimensional (QID) orbital fluctuations [16]. 

2. Magnetic and orbital order in cuprates and manganites 

Conceptually the simplest realistic spin-orbital model can be derived for 
ions interacting on a cubic lattice, as in KCUF3. The charge excitations 
dfd^ # d^dY lead to: one high spin state, and two low-spin "^E and 
^Ai states [8]. The energy spectrum in Fig. 1(a) is obtained from the 
model Hamiltonian which includes the on-site U and Jh interactions for 
degenerate d orbitals [17], and reproduces the exact spectrum [18]. The 
superexchange is oc Jg = t^/U, where t^ is the largest hopping element 
between two — orbitals along the c axis (note that this is a natural 
unit for the anisotropic hopping between eg orbitals [5]), and is given by 

n{d^) = JY: E [{s.-s, + s^)ji]\d') + k\]\d% (1) 

where Si arc spin S = 1/2 operators. The operator expressions: 

Jlf (d^) = (2 + np2 - r?P3)^g) - riiSpi - p2)Pg.) , (2) 
k^]\d') = -[1 + r?(3pi +P2)/2]Vg - [l + v{P2-P3)/2]Vg, (3) 

describe spins and orbital superexchange, with 77 = Jh/U, pi = 1/(1 — 3r]), 
P2 = 1/(1 — r]), and pa = 1/(1 -|- 77). They depend on orbital operators: 

Vg = (1/2 + r7)(l/2 - r]) + (1/2 - r7)(l/2 + r/), (4) 
Vg = 2(1/2 -r7)(l/2-r7), (5) 

which project on the orbital states, being either parallel to the bond (ij) 
direction on one site (Pj^ = 1/2 — t^) and perpendicular on the other {Pj^ = 
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1/2 + rj), or parallel on both sites. They are represented by the orbital 
operators associated with the three cubic axes (7 = a, 6, or c), 

= {-at ± V^af)/^, rf = at/2, (6) 
where the ex's are Pauli matrices acting on: jx) = ( 3. ), |z) ^ 



^ ' ' VI 
which transform as |a;) oc a;^ — 

The superexchange in LaMnOs couples total spins 5 = 2 at the Mn^+ 
ions and originates from the charge excitations, dfd^ # dfd^ [9]. The eg 
part, following from ^ (^2g)'^j (^Ig^g) processes, involves FM terms 
due to the high-spin ^Ai state, and AF terms due to the low-spin states: 
'^^i, '^E, and [Fig. 1(a)], is orbital dependent. By contrast, the t2g 
part oc it ~ 0.09, which follows from dfd^j ^ (^2g)^j(*2p^ff) excitations, is 
purely AF and orbital independent. Both terms give 

nd') = j,Y. E [{Si-s,+^)j\]\d')+k\]\d% (7) 

where the exchange interactions depend on the multiplet structure. 



= 4l-t^(^3 + 2g4)]7'« -^7?(36gi+9g2 + 20g3)7'g)+jt, (8) 
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with qi = 1/(1 - 3r/), 92 = 1/(1 + 2r?), 93 - 1/(1 + 8r//3), and 94 ^ 
1/(1 167?/3) [18]. The orbital part K^f {d^) is given in Ref. [9]. 

Both d? model and model at jt = describe strongly frustrated 
superexchange in the limit of Jh 0, which takes a universal form, 

=JeY.T. p. ■ Sj/S' + l)(l/2 - r7)(l/2 - rj) - l] . (9) 

Several classical phases have the same energy of —3 J,, per site at this point 
[6]: the G-AF phases with arbitrary occupation of orbitals, and A-AF phases 
with ((1/2 — t7)(1/2 — rJ)) =0, as obtained for staggerred planar orbitals, 
e.g. for — ip' jy^ — z^. The model (9) is qualitatively different from the 
idealized SU(4)-symmetric case [19] due to the directionality of eg orbitals. 

At finite Jh the degeneracy of classical phases is removed, and the ^-AF 
phase is stable, with two-sublattice alternating orbital order in both cuprate 
(1) and manganite (7) model, = cos 9 i\i z a) ^svn 9 i\ix a), where it refers 

to i G A{B) sublattice. In the cuprates the orbital order given by cos 2^ = 
(1 — 77/2)/(2-|-377), induces FM interactions Jah within the (a, 6) planes, and 
AF interactions Jc between them [5]. The AF interactions decrease with 
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increasing Jh/U [Fig. 2(a)], but still dominate at realistic Jr/U ~ 0.12 
[7], explaining why the excitation spectra of KCUF3 are dominated by QID 
spin excitations of S = 1/2 spin chains [20]. 

Although the orbital order found in the manganitc model (7) at Jh/U = 
is again — /y^ — , and the A-AF phase is stable, the situation is 
here qualitatively different as Jab and Jc change much faster with increasing 
Jh/U [Fig. 2(b)], and have similar values in LaMnOs {Jh/U ~ 0.117 [7]), 
demonstrating the proximity to ferromagnetism which is indeed observed 
in doped manganites [1, 5]. Including the (smaller) t2g interactions one 
finds a somewhat enhanced tendency towards antifcrromagnetism, with the 
G-AF {A-AY) phase stable for Jh/U < 0.05 {Jh/U > 0.05). In order to 
explain quantitatively the experimental ratio Jc/ Jab — 0-7 in LaMnOa, one 
has to include also the JT effect which stabilizes the orbital order closer to 
(|a:) + |z))/(|a:) — |z)) alternation [9] . This modification of the orbital ordering 
changes not only the effective magnetic interactions, but also considerably 
reduces the scattering of a hole on spin excitations in LaMnOs [21]. 



3. Orbital fluctuations in t2g systems 

As in the case, the excitation spectra of and ions in the t2g 
subspace [18], shown in Fig. 1(b), may be faithfully reproduced with a 
model Hamiltonian [17] containing only two parameters: U and Jh, with Jh 
standing now for the Hund's element between two t2g orbitals. As usually, 
the excitation energy to high-spin (^^2 and ^^2) states \s U — 3Jh, while 
the energy of the next (low-spin) excited states is either U — Jh for d^ ions 
(^T2, ^E), or U for d^ ions (^Ti, '^E), respectively. The highest excitation 
energy of C/ + 2Jh is the same for c?^ (""^^i) and d^ CT2) ions. 

Each t2g orbital is orthogonal to one of the cubic axes, so we label them as 
a, b , and c (for instance, xy orbitals are labelled as c). The superexchange 
interactions oc J = /U follow from the hopping between two orbitals 
active along a given direction 7, for instance between the pairs of a and h 
orbitals along the c axis. Therefore, it is convenient to define pseudospin 
operators, % = {Tf,Tf, rf }, which act in the subspace spanned by two active 
orbital flavors [10, 16]. For instance, for a bond (ij) \\ c, these operators 

are: r+ = af^^, r" = 6ja^, rf = ^{n^^ - n^^), and nP = n^^ + n^^, where 
{aj,6|} are Schwinger bosons for a and h orbitals. 

The model for titanates follows from the d\dh ^ ti^d^ processes. 



n{d') = jY. E {{Sr-s,^s')jtl'{d^)^kt^^{d% (10) 
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with the exchange constants between S = 1/2 spins, 

jI'^^ = 2{Ti-Tj-ininj/4)+ri[{-3ri+r2){nianjb+nibnja+nic+njc^ 

+ (3ri+r2)(r+r-+r-r/) + 4(r2-r3)(nian^„+ Wj6)/3]/2, (11) 

depending on: ri = 1 — 3r/, r2 = 1 — ry, ra = 1 + 2rj, while K^j\d^) 
stands for purely orbital interactions. A priori, the magnetic interactions 
are anisotropic, and may be cither AF or FM, depending on the orbital 
correlations. In the limit of Jr/U = the Hamiltonian (10) takes the form, 

= (J/2) ^ [^Si-S^/S^ + l){n-rj + n^n^l^) - 45/3], (12) 

and shows again a strong frustration of superexchange interactions [10]. 
Although formally it resembles the SU(4)-symmctric spin-orbital models [19] 
even more than Eq. (9), the pseudospin operators Tj have here a different 
meaning and refer to different orbital flavors for each cubic direction 7. 
One may also notice a certain analogy with the models of valence bond 
solids [22] , but this analogy is again only partial, as the formation of orbital 
singlets in all directions simultaneously is impossible. 

In the mean field approach (MFA) the G-AF phase is degenerate with 
FM phases, if (rj-'r, + jriirij) = 0, as realized for alternating orbitals (e.g. 
for staggered a/b orbitals). Such FM states, with anisotropic exchange 
constants: Jpa and Jpc along a (b) and c axis [Fig. 3(a)], respectively, 
would be favored classically at finite Jh- On the contrary, the quantum 
fluctuations take over, remove the anisotropy, and stabilize the orbital liquid 
state, if the JT interactions are weak [10]. Indeed, the spin wave spectrum 
of LaTiOs is nearly isotropic [11], showing that the orbital moments of t2g 
ions are fully quanched [10]. Increasing Jh almost does not change the 
exchange constants Jaf evaluated using the MFA in this state [Fig. 3(a)]. 

The superexchange interactions between S = 1 spins in LaVOs [16], 

n{d') = jj: E [iS^■S, + l)J^J\d') + K^]\d% (13) 

T fe>ll7 



follow from the dfdj ^ djd^ processes active on the bonds, with 



{l + 2r]R)(Ti-fj+^nin^ -rjr^rlr^+^nin^ - ^r]R{ni+nj) 



(7) 



(14) 



and the orbital term K^^'j'' given in Ref. [16]. The coefficients R = 1/(1 — 37/) 



and r = 1/(1 + 27/) follow from the multiplet structure of d^ ions [Fig. 1(b)]. 
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In the limit of J// — > one finds again the frustrated superexchange (12). 
While the orbital liquid cannot stabilize in this case, orbital singlets may 
form along the c direction when c orbitals have condensed (n^c = 1) and the 
o and b orbitals fluctuate. This gives a novel mechanism of ferromagnetic 
interactions which operates already in the limit of Jh = [16]. 

The exchange constants within (a, 6) planes (Jab) and along c axis (Jc): 

Jab = [l-viR+r) + {l + 2r]R-rjr){nian,a)^''^]/4, (15) 
Jc = [{l+2rjR){n-rj + l/4)(^) - 77r(r/r/ + l/4)(^) - nR]/2, (16) 

are given by orbital correlations. Their values at r/ = were obtained 
from the Bethe ansatz for a QID Heisenberg chain, while the orbital wave 
spectrum, = [A^ + R^{1 - cos^ k)]'^^^, with a gap A = {r]{R + r)[2R + 
r]{R + r)]}^/^, was used at finite Jh- As a result, one finds increasing FM 
(Jc) and decreasing AF (Jab) exchange constants with increasing Jfj [Fig. 
3(b)], and both interactions have similar values at Jh/U ~ 0.15 [7]. 

While the cubic structure of LaVOs is almost undistorted [14], YVO3 
has a distorted structure, and a and h orbitals stagger in (a, 6) planes and 
repeated themselves along c axis [15]. Such ordering can be promoted by 
the JT effect term which lowers the energy by —2V on the bonds along 
the c axis when a (b) orbitals are repeated in C-type orbital ordered state 
[16]. Finite V > lowers the energy of the G-phase, but the entropy S 
determined by orbital excitations increases faster in the C-phase, and thus 
induces a transition from G-AF to C-AF order around T* ~ 0.8 J (Fig. 4), 
reproducing qualitatively the first order transition observed in YVO3 [15]. 

4. Summary and open problems 

In summary, the transition metal oxides with orbital degrees of free- 
dom show a very fascinating behavior, with various types of magnetic and 
orbital order. While Cg orbitals usually order and explain A-AF phases, 
further stabilized by the JT effect, the t2g orbitals have a generic tendency 
towards disorder, which leads to the orbital liquid in the isotropic G-AF 
phase in LaTiOs. In cubic vanadates the JT interactions compete with the 
orbital disorder, and the QID orbital fluctuations stabilize the C-AF phase 
in LaVOa, and also in YVO3 at finite temperatures. A better understand- 
ing of these fluctuations is required to explain quantitatively the observed 
phase transitions and the strong reduction of the magnetic order param- 
eter in LaVOa and YVO3. This problem is as urgent as the theoretical 
understanding of the colossal magnetoresistance in the manganites. 
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Fig. 1. Excitation spectra in cubic transition metal oxides for: (a) eg systems: 
Cv?+ {d^) and Mn'^+ {d^) ions; (b) t2g systems: Ti'^+ {d^) and (d^) ions. 



Fig. 2. Exchange constants FM Jab (solid lines) and AF Jc (dashed lines) in A- 
AF phase of Cg systems as functions of Jh/U for: (a) cupratcs (KCuFs); (b) 
manganites (LaMnOs), for: jt = (thin lines) and jt = 0.09 (heavy lines). 



Fig. 3. Exchange constants as functions of Jr/U for t2g systems: (a) G-AF {Jaf, 
dashed line) and FM {Jpa and Jpc, solid lines) phase in titanates; (b) AF J^b 
(dashed line) and FM Jc (solid line) C-AF phase in vanadates (LaVOa). 



Fig. 4. Free energies F{T) ^ {n{d^)) - TS (in units of J) of: G-AF phase obtained 
with the JT interaction V = 0.65 J (sohd line), and G-AF phase for ry = 0.05, 0.10 
and 0.15 (dashed lines), as functions of temperature T/J (after Ref. [16]). 



U+5Jh^ 



(a) 4^, 



U+2J 



H 



- 1 



A 



U-JhI- 



U-3J H I- 



•A 



|8 



U+2J H I- 



U-JhI- 
U-3J H I- 







'2 






T2^E 











